Light scattering from ultracold atomic gases in optical lattices at finite temperature 
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We study light scattering from atoms in optical lattices at finite temperature. We examine the 
light scattered by fermions in the noninteracting regime and by bosons in the superfluid and Mott 
insulating regimes. We extend previous theoretical studies to include the full band structure of the 
optical lattice. We find that light scattering that excites atoms out of the lowest band leads to 
an increase in light scattering away from the classical diffraction peaks and is largely temperature 
independent. This additional light scattering leads to lower efficiency of temperature measurements 
based on photon counting. 
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Ultracold atoms in optical lattices exhibit a remark- 
able range of behavior and have become an indispensable 
tool for investigating many-body phenomena [HQ- As 
the gamut of many-body states achieved in experiments 
with optical lattices expands, it is also important to con- 
sider how to detect and reveal more information about 
these states. With this goal in mind, many authors have 
considered the scattering of light by these ultracold gases 
as a potential way of extracting information. 

Bragg scattering, where light scattering from the sam- 
ple is stimulated by a pair of lasers, is a prominent 
method that has been investigated as a probe of bosons 
in optical lattices 3-7], and experiments using Bragg 
scattering have detected signatures of the superfluid and 
Mott insulator states of the Bose-Hubbard model || Q . 
Another suggested scheme involves light scattering into 
a cavity, where the strong light-matter coupling encodes 
information about the many-body state into the distri- 
bution of photons transmitted by the cavity flol - fTij . The 
effect of finite temperature on light scattering has been 
investigated as a potential probe of temperature for the 
superfluid and Mott cases [15| and also for the noninter- 
acting Fermi-Dirac gas [l6| . Recently the far field diffrac- 
tion pattern of a Mott insulator was imaged in an experi- 
ment where the diffraction peaks were used to detect the 
spin distribution fl7j . Detection of spin distributions in 
optical lattices using light-matter interactions has also 
been considered theoretically (l8l - l27| . 

In this paper we consider what light scattering, which 
is not stimulated by another laser or by a cavity, tells 
us about the many-body state. The aim being to re- 
veal information about the atoms while they remain in 
situ with a relatively simple experimental implementa- 
tion, only requiring a laser to be shone on the atoms, as 
depicted in Fig. [TJ followed by the detection of scattered 
photons, possible after collection by a lens. This would 
remove the need for a cavity or the need to destructively 
image the atomic sample, which is usually required with 
Bragg scattering (see Ref. [28| for a recent exception). 
When scattering is not stimulated there is less restric- 
tion on the final energy state of the scattered photons 
and the full band structure of the optical lattice must be 
taken into account, something which has been neglected 



FIG. 1. (Color online) Light scattering scheme. The atoms 
in the optical lattice are illuminated by a laser beam with 
wave vector ki. Light is then scattered to other modes k/ 
with spherical coordinates 9 and <j>- In our examples of light 
scattering the optical lattices are two-dimensional square lat- 
tices with lattice translation vectors ax and ay, and the input 
light has ki = nz/a. In this coordinate system we then have 
k = ki — kf — — 7r(sin#cos</>, sin sin 0, cos 9 — l)/a. 



in previous works. We present here a treatment of light 
scattering from ultracold atoms in optical lattices that in- 
cludes both the band structure and the effects of nonzero 
temperature. 

Our approach is organized as follows. We begin by 
describing the light scattering process that forms the ba- 
sis of our work. We then examine the light scattering 
patterns that result for noninteracting fermions, and for 
bosons in the superfluid and Mott insulator regimes of the 
Bose-Hubbard model. We then examine how the temper- 
ature measurement described by Ruostekoski et al. [lfjj 
applies to these systems and is affected by the multiband 
structure of the optical lattice. 
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I. FIRST ORDER LIGHT SCATTERING 



mode (k/,A) is then 



Nonrcsonant light scattering from an ultracold gas can 
occur in two ways. The first involves the diffraction of 
light due to the nonuniform density of the gas, which 
leaves the gas unchanged, while the second involves the 
creation of excitations in the gas where the light transfers 
momentum and/or energy to the gas. In the following wc 
see how these processes occur when laser light is scattered 
from atoms in first-order perturbation theory. 

We consider light-matter interactions where the atom 
can be treated as a two-level system, as may be the 
case because of angular momentum considerations in the 
dipole approximation. For example, in 87 Rb, the transi- 
tion \5 2 S 1/2 ,F = 2,m F = 2) -> |5 2 P 3/2 , F = 3, m F = 3) 
that is excited by cr + -polarized light is a cycling transi- 
tion, where, because of the polarization of the input light, 
only one excited state can take part, which in turn can 
only decay to the initial ground state. Further simplifi- 
cation can be made when the light incident on the atoms 
is off- resonant, with detuning A = u a — u>i between the 
atomic resonant frequency u> a and the input laser fre- 
quency uji. In which case the excited atomic state can 
be adiabatically eliminated, giving, in the rotating wave 
approximation, the following interaction Hamiltonian, 



H 



dr¥(r)E~(r) ■ E+(r)£(r), 



(1) 



where g = ineoc 3 ^ /(ui^A) and 7 is the lifetime of the 
excited state. The Hamiltonian is expressed in terms of 
the ground-state atomic- field operators ^(r) and &(r) 
and the positive and negative frequency components of 
the electric field, E + (r) and E~(r). 

In the interaction picture we have Hj{t) = 
e iH t/hjj e -tH t/h^ wnere jj i s the Hamiltonian of the 
noninteracting light-matter system, and the time evolu- 
tion operator is given by the integral equation Ui(t, 0) = 
1 - |/o dtfHiWU^tf, 0) HI]. For a weak perturbation 
the evolution of the state is given approximately by ex- 
panding this integral equation to first order, giving 



!*(*)>; 



dt'Hj |*(0))j 



(2) 



We take the initial state of the system to be an atomic 
cigenstate \^f u ) and a laser represented by the classical 
field E + (r) = |e lki ' r e.;. The first-order coupling then 
leads to photons being scattered from the laser into other 
modes (k/,A) with wave vector kf, frequency u)f and 
polarization e,\(k/), resulting in a momentum change 
fik = h(k{ — kf) and energy change hw = h(u>i — u>f). 
The probability at time t of the system being in a new 
atomic eigenstate \*& v ) with an additional photon in the 



k/W 



^ Jo 



dt'H! 



xS t ((E v -E u )/h-Lo) 



dre ik - r (*„|^(r)*(r)|* u ; 



(3) 

where G x {k f ) = \g£e* x {k f ) • e 4 | 2 W/ /(8fte (2^) 2 ) is the 
coupling constant between the electromagnetic field 
modes and Ej is the unperturbed energy of the state 
The function S t (oj) = 2sin 2 (wi/2)/(7rtaj 2 ) approaches the 
Dirac delta function 5(oj) as the interaction time t ap- 
proaches infinity, enforcing energy conservation. 

For cr + -polarized input laser light we can sum over po- 
larizations in the coupling to get G(k/) = J2\ ^(k/) = 
\g£ 1 2 ( 1 + cos (60 2 )u)f I ( 16he (2tt) 2 ) . The total rate of scat- 
tering photons with wave vector kj and frequency uif of 
cither polarization is then 

r(k />W/ ) = G(k / )5(k J - k /jWi - u>f), (4) 

where S(k,u) is the dynamic structure factor [13] ■ For 
a finite-temperature system, in the canonical ensemble 
wi th p artition function Z, the dynamic structure factor 

is MM 



s(k,u) = jJ2 6 « e « ^y~ Eul{kBT 



dre ik ' r (^|*t( r )|r( r )|^ 



(5) 



The structure factor can be divided up into two parts, the 
first with u = v describes classical diffraction, which re- 
sults in no energy or momentum transfer to the gas, that 
is, the light scattering is elastic. The second part, where 
h/d. describes inelastic light scattering that results in 
excitations of the gas. 

For nonresonant light scattering the frequency change 
lj is determined by the difference in energies between the 
initial and final atom many-body states. For scatter- 
ing from ultracold atoms these energy differences are all 
many orders of magnitude less than the frequency of the 
input light. To a very good approximation we then have 
kj = |kj|k/, where kf has associated angles 8 and as 
shown in Fig. [TJ To get the total rate of photon scatter- 
ing in direction kf, we integrate Eq. Q over frequency. 
This gives 

Itotalfcf) = Iatom(0) X S(k, - \k.\kf), (6) 

where I atom (8) = 9/ m7 2 (l + cos(9) 2 )/(32hckf A 2 ) is 
the scattering distribution resulting from the electronic 
structure of the atom and laser intensity Ij n = eoc|£| 2 /2, 
and 

S(k) = \Y.\ drdv'e^-^e-^/^ 

U 

X (* u |*t( r ')^( r ')^rt( r )4( r )|^ u ) (7) 
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is the scattering distribution due to the spatial structure 
of the atomic sample, known as the static structure factor 

The static structure factor applies in the perturbative 
regime, where it has been used successfully to describe 
neutron scattering from liquid helium [3(| and to calcu- 
late the spectrum of light scattered from ultracold gases 
[HI, HU, along with the response of ultracold gases to 
Bragg excitation [34[. In this paper we use the struc- 
ture factor to predict the angular dependence of light 
scattering for atoms in an optical lattice. While our for- 
malism can be applied to any lattice and input laser con- 
figuration, in all our examples we consider the setup in 
Fig. [U where the optical lattice is assumed to be a two- 
dimensional square lattice in the xy plane, illuminated 
by light propagating in the z direction with wavelength 
equal to that of the light used to create the optical lat- 
tice. The change in photon wavevector upon scattering is 
then k = — k/ = — 7r(sin# cos </>, sin 9 sin0, cos 9 — 1) /a, 
where <j> and 9 are the scattering angles shown in Fig. Q] 
and a is the lattice intersite separation. We denote the 
strength of the optical lattice potential in the x, y and 
z directions by V x , V y , and V z , respectively, in units of 

the recoil energy Er — ^ ^ 2 • ^ n ^ ne following sections 
we calculate the light scattering predicted for fermions 
and bosons and see how different correlations in different 
many-body states affect the angular distribution of the 
scattered light. 



II. FERMIONS 

For spin-polarized fermions of the same species, s-wave 
scattering is prohibited, and for low temperatures the 
fcrmionic atoms are approximately noninteracting. The 
Hamiltonian for the many-body system is then exactl y di - 
agonalized by the single-particle Bloch states ^ q , m (r ) |35| , 
where q is the quasimomcntum and m is the band in- 
dex. At finite temperature each Bloch state is occupied 
according to Fermi-Dirac statistics. 

We can expand the atomic-field operators in terms of 



the Bloch states, *(r) = m b 



q.mrq.m 



(r), where 6 qv 



is the annihilation operator of an atom in the Bloch state 
0q,m(r). Together with the corresponding set of cre- 
ation operators 6 q m , these operators obey the standard 
fermion anticommutation relations. This expansion leads 
to the following expression for the static structure factor, 



,mi,q 2 ,ni2 



q4,m4,q3,m3 



(k) 



qi,q2,q3,q4 

mi ,rxi2 ,1113 ,iri4 



(^qi,mi^q2,m 2 ^q 3 ,m3^ 



q3,lT13 "q4,m4/ ; 



(8) 



where 



/qi,m 1 ,q 2 ,m 2 

(k) = M J <ir<£* 2 m2 (r) 



^(r)e l * r (9) 



is the transition matrix element between the Bloch func- 
tions associated with a momentum transfer k. For suf- 



ficiently low temperatures all the atoms reside initially 
in the lowest band of the optical lattice; however atoms 
may be excited to higher bands by scattering a photon. 
We then have two terms, the first due to light scattering 
within the lowest band of the optical lattice, that is, 



S 9 (k) 



qi,q2,q3,q4 



/qi,0,q 2 ,o(k)/q 4 ,0,q 3 ,o(k) 



X (&q 1 ,0&q 2 ,0&]j3,0 b q4,0), (10) 



and the second due to scattering into the higher bands 
of the lattice, that is 

S&(k) = /qi,0,q 2 ,m(k)/q 3 ,0,q 2 ,m(k)(S qii0 S q3 ,o). 

qi,q2,q3,m7io 

(11) 

This part of the structure factor was neglected in other 
treatments of finite-temperature light scattering [ill, |T(| , 
and we find that this can make a significant difference to 
predictions about the light scattering. 

We can calculate the higher band component with the 
help of the sum rule 



•/'qi,0,q,m(k)/q2,0, q ,m(k) — S, 



qi,q2i 



(12) 



q,m 



which is implied by the completeness of the Bloch func- 
tions. Using this relation we find 

S b (k) = N- £ /q 1 ,0,q 2 ,o(k)/q 3 ,0,q 2 ,o(k)(S qi)0 S q3j o), 

qi,q2,q3 

(13) 

where N is the total number of atoms on our M-sitc 
lattice. All operators are now in the lowest band and we 
drop the band index from our notation. 

We can now evaluate all the expectation values us- 
ing Fermi-Dirac statistics. We find that light scattering 
within the lowest band is the sum of two components 
<S 9 (k) = S g o(k) + SgiQs). The zeroth-order component, 
which does not result in many-body excitations in the 
lattice, is 



S s o(k) 



-£|iq,q(k)| 2 7V q (l- 



(14) 



where 2V q = (6 q 6 q ) gives the expected population of the 
lowest band Bloch state with quasimomcntum q. This 
component is the sum of a contribution due to the av- 
erage density of the Fermi-Dirac gas and a contribution 
due to the finite-temperature fluctuations in the average 
density. The first-order component, resulting from atoms 
being excited into unpopulated Bloch states within the 
lowest band, is 



S gl (k) = ]T |/q,p(k)| 2 7V q (l - N p ). (15) 

q^p 
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FIG. 2. (Color online) Comparison of the angular dependence 
of the components of the static structure factor for a Fermi- 
Dirac gas on a 30 x 30 lattice with V x = V y = 8Er, V z = 15Er 
and temperature 0.02_Er/£;b = 0.16Tf, where Tf is the Fermi 
temperature, (a) Dependence on 6 of the structure factor 
components 5' 9 o(k) (solid blue line), S 9 i(k) (dashed green 
line) and St (k) (dotted red line) for cj> — and filling factor 
/ = 0.5. (b) As in panel (a) with <f> — tt/4. (c) Comparison 
of the structure factors for gases with filling factors / = 0.5 
(solid blue line) and / — 1 (dashed green line). 



The scattering function for scattering involving the 
higher bands, which we refer to as interband scattering, 
is given by 



S fe (k) = iV-£|/ q , p (k)| 2 iV q 



(16) 



The above theory extends the treatment in Ref. [16[ 
through the inclusion of the interband scattering and, 
further, by consistently working in the Bloch basis we 
have not made any approximations involving the local- 
ization of atoms to lattice sites, allowing our model to be 
applied over the full range of lattice strengths. 

The components of the light scattering are shown in 
Fig. [2] as a function of scattering angle for a two dimen- 
sional lattice as described in Fig.[TJ The strongest compo- 
nent is the elastic scattering, which forms sharp diffrac- 
tion peaks with a diffraction maximum at 8 = 4> = 



that scales with N 2 . The interband component scales 
with N and is near zero intensity at the diffraction maxi- 
mum, but can overwhelm the elastic scattering for higher 
angles and leads to a background level of scattering be- 
tween the classical diffraction peaks. Neither the elastic 
term nor the interband term are significantly affected by 
the temperature or occupation statistics of the gas. The 
first-order component does change however, for example 
if the lowest band is completely full with N = M, no 
excitations are available within the lowest band due to 
Pauli blocking and the first-order component is identi- 
cally zero, while for N/M = 0.5 the contribution from 
this component scales with N. In Fig. [2jc) we compare 
the angular distribution of light scattered for a full lattice 
with that of a lattice with filling factor / = N/M = 0.5. 

The first-order component also decreases as 9 — > due 
to the Pauli blocking that appears in Eq. (fTBf . Fermions 
can only scatter into states that are unoccupied and for a 
particular momentum transfer Kk this restricts the num- 
ber of atoms that can scatter light. At zero temperature, 
the only atoms that can be involved in inelastic scatter- 
ing within the first band arc those with momenta hq such 
that ft|q + k| > pp, where pf is the Fermi momentum. A 
lower scattering angle corresponds to a lower momentum 
transfer, reducing the number of atoms that can scat- 
ter light and hence the amplitude of light scattering. As 
temperature increases more final momentum states be- 
come available to scatter to throughout the momentum 
distribution and the low-angle scattering increases with 
temperature. We further examine the temperature de- 
pendence of the light scattering in Sec. IIV1 



III. BOSONS 

For bosons in an optical lattice at low temperatures s- 
wave interactions between the atoms play an important 
role. The resulting system d yna mics are well described 
by the Bose-Hubbard model [36j], which is characterized 
by the energy of on-site interactions U and the energy 
associated with tunneling between lattice sites, J. For 
low on-site interaction strength, the atoms form a super- 
fluid. As interaction strength increases a phase transi- 
tion occurs and the system becomes a Mott insulator. 
We examine the light scattering from both these phases, 
beginning with the superfluid regime. 



A. Superfluid 

In the superfluid regime the atoms are not localized 
at lattice sites and, particularly for light scattering, it 
is most intuitive to expand the atomic field operators in 
the Bloch basis as wc did for the noninteracting fermions. 
The creation and annihilation operators for the Bloch 
states now obey the standard boson commutator rela- 



tions. The atomic Hamiltonian is then 

h = Y, e ^ + \ E 



f^qi .q2 ,q3 ,q4 k qi ^q 2 ^q3 ^q4 j 



qi,q2,q3,q4 



(17) 



where 



17. 



4irh 2 acAI 



qi,q2,q3,q4 



rfr^q! ( r )</>q 2 ( r )^q 3 ( r )0q4 0), 

(18) 



results from the s-wave interaction with scattering length 
a s . Here for notational convenience we have assumed the 
band index of the Bloch states is contained in the gener- 
alized quasimomentum q, and the theory below includes 
all bands. 

At zero temperature, atoms in a non-interacting Bosc 
gas condense into the zero momentum Bloch state 4>o(r). 
For weak interactions this ground state is perturbed, and 
the new ground state and elementary excitations of the 
gas are found using Bogoliubov's theory [37J. Extensions 
of Bogoliubov's theory to Bose gases in optical lattices 
have been made [38l |39| , but work in the tight binding 
approximation. Here, by working in the Bloch basis we 
do not need to make this approximation. 

Following Bogoliubov's treatment we replace the oper- 
ators bo and b of the highly populated zero momentum 

(b Q b ). Then as- 



Bloch state by the c-number t/INo 
suming the nonzero-momentum Bloch states each have a 
population much smaller than No, we can neglect terms 
above quadratic order in the Hamiltonian. The Hamilto- 
nian can then be diagonalized via the Bogoliubov trans- 
formation (37| 

= U q/3q - «q/3-q, 6 q = M q /3 q - Uq/3I q , (19) 

where the Bogoliubov transformation coefficients are 
given by 



\ 




and 



\ 



2 \ too 



(20) 



with Eq = £ q - E - N Q Uo + 2AW q and C/ q = C/q,o.o,q- 

The operators /3 q and /3 q are interpreted as the creation 
and annihilation operators of noninteracting bosonic 

quasiparticles with energy hujq - 



iV 2 [/2. At zero 



temperature the ground state of the system is the quasi- 
particlc vacuum, while for finite temperature T the pop- 
ulation of the quasiparticle modes is given by the Bosc- 
Einstein distribution [HI H~Dl 



q.q' 



= 8, 



q,q' n q- 



(21) 



Furthermore, we have (/? q ) = 0t) = and 

(/? q /3 q /3 q ,^q') = n q n q / + £ q , q '(n q + l)n q , which we will 
need to calculate the structure factor. 
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FIG. 3. (Color online) Comparison of the angular dependence 
of the components of the static structure factor in the super- 
fluid regime for a 30 x 30 lattice with 2700 atoms. The lattice 
strengths are V x = V y = SEr and V z — 20Er, leading to 
parameters J = 0.11i?R and U = 0.17Er, the temperature is 
0.05i?fl/fcfl and the condensate population is No = 2527. (a) 
Dependence on 9 of the structure factor components S g o(k.) 
(solid blue line), S 9 i(k) (dashed green line), Sg^Qs) (dot- 
dashed black line) and S&(k) (dotted red line) for <f> = 0. 
(b) As in panel (a) with <j> — n/4. (c) Inelastic components 
only. 



The above quantities depend implicitly on the conden- 
sate number No, which is in turn restricted by the fol- 
lowing relation for the number operator N = ^ c btb G 



jq ,y q u q' 



(7V)-^iVq = 7Vo 



E 

q#0 



u\n c 



1)) 



(22) 



For a fixed (N) = N, Eqs. (JSOj) and ([22]) must be solved 
self-consistcntly. 

As noted above, the Bogoliubov approximation only 
includes the contributions to the interaction energy pro- 
portional to Nq, and its validity is therefore limited to 
the regime where the number of quasiparticles is a small 
fraction of No- In this work we limit our calculations so 
that this fraction is less than one tenth. 

As in the fermion case, we assume the temperature is 
sufficiently low so that all the atoms initially reside in 
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the lowest band of the optical lattice; however atoms 
may be excited to higher bands by scattering a pho- 
ton. Examining Eq. (|20|) we see that if, for the higher 
bands, Eq 3> Ngllq, then u q ~ 1 and the quasipar- 
ticle excitations arc essentially particle-like. This cri- 
terion is easily satisfied in the cases we investigate, in 
which case the only nonzero averages in the static struc- 



^1,0^2,0^3,0^4,0; 



and 



ture factor given by Eq. (JSj) are (b 

(& qi ,0^q2,m& q2 , m f>q 3 ,o) = (& qi , &q 3 ,o) for m ^ 0. 

We can now evaluate all the expectation values us- 
ing the Bogoliubov theory for the lowest band, and 
the quasimomenta q from this point refer only to the 
lowest band. We find that light scattering within the 
lowest band is the sum of three components S g (k) = 
5 s o(k) + S' 9 i(k) + 5 92 (k). The zeroth-order component, 
which does not result in excitations in the lattice, is 



Ew k )^ 

q 

+ E l/q,q( k )l 2 K + «q) 2 "qK + !)• (23) 



As in the fermion case this component results from the 
average density of the gas and the finitc-tcmpcraturc fluc- 
tuations in the average density. The first-order compo- 
nent, resulting from atoms being excited out of or into 
the condensate mode, is 

S gl (k) =N J2 l/o, q (k)| 2 («, - vqf (2n q + 1), (24) 
q/o 

and the second-order component is 

S 92 (k) = E l/q,p( k )| 2 ( w q"p + VqV p ) 2 riq(n p + 1) 
q#o, P #o 
q#p 

+ \ E l/q,p( k )| 2 ( u q v p +Vau p ) 2 
q#o, P ^o 

x (n q n p + (n q + 1) (n p + 1)) , (25) 

which results from creation and/or destruction involving 
two quasiparticlc modes. The distribution resulting from 
interband scattering is given by 



S h (k)=A-^|/ q , p (k)| 2 7V q . 



(26) 



Again, our theory extends that of previous studies |6lll5l| 
by including the interband scattering and by working 
consistently in the Bloch basis, with no localization ap- 
proximations, we have a model that is applicable even 
for very weak optical lattices. 

In Fig. [3] we compare how the four components of the 
static structure factor vary as a function of scattering 
angle. We see that the classical diffraction pattern is 
again the dominant feature with peaks of order N 2 , while 



the first order and interband components are of order N. 
The second-order component is larger than the first-order 
and interband components near 8 = 0, a behavior that 
persists even at T = 0. 

The behavior of the first- and second-order components 
reflects the correlations between atoms in the lattice and 
quantum interference and enhancement play a role in de- 
termining their structure. We sec in Fig.[3]that the first- 
order component vanishes as k — > 0. This phenomenon 
has been recorded in Bose liquids [HI, E2 and in an ex- 
periment with a Bose condensed dilute gas, where Bragg 
spectroscopy was used to probe the structure factor [43| . 
Here the same processes are at work and we can un- 
derstand the vanishing of the first-order component as 
interference between two scattering channels. The first 
channel involves the quantum depletion of the conden- 
sate, which consists of pairs of atoms with opposite mo- 
mentum for finite intcrparticle interactions [44|. Light 
scattering can scatter an atom in one of these pairs back 
into the condensate. This gives the same final state as 
created by the second scattering channel where an atom 
is scattered out of the condensate mode, and the ampli- 
tudes of these scattering processes destructively interfere 
as k — >• (see Ref. [43[ for further details). An interest- 
ing feature of the optical lattice case, as opposed to the 
uniform Bose gas, is that the rate at which the first-order 
structure factor goes to zero depends on the ratio of J to 
U, which can be tuned by adjusting the height of the lat- 
tice or by Feshbach resonance. By adjusting the lattice 
parameters and examining the light scattering at various 
angles we can see how the depiction, or alternatively how 
number squeezing [3^ |. in the lattice changes. 

Despite involving two quasiparticlc modes, the second- 
order component can still make a significant contribution 
to the light scattering. This is due to quantum enhance- 
ment in the higher order processes as seen in the ex- 
periment by Rowen et al. [45j . For this component to be 
nonzero, noncondensate modes must be populated, either 
as a result of intcrparticle interactions or finite tempera- 
ture. 



Mott insulator 



Wc now examine light scattering from the Mott insula- 
tor. This state is characterized by strong interactions be- 
tween atoms leading to particle localization and one may 
expect that the characteristics of light scattered from this 
state will be significantly different to the superfluid case. 

Again we assume that the temperature is low enough 
that all the atoms initially reside within the lowest 
band of the lattice. Expanding the atomic-field opera- 
tor in terms of the localized Wannicr site basis, ^(r) = 
Y]. m bj^ m w m (r — Yj), gives the following expression for 



the static structure factor 

,0,33 ,m (k) 

jl J2 J3,j4,m 

X ''i - • (27) 

where 

/ji,mj a ,n(k) = /" dr Wm (r - r^)u; n (r - r j2 )e lk r . (28) 

The overlap between lowest-band Wannicr functions at 
different sites is very small in the Mott regime and 
/n,o,j2,o(k) can be neglected for j\ fa- 

As in the previous cases the terms involving higher 
bands can be simplified using a sum rule, this time for the 
Wannicr functions. Once again, the description now only 
involves contributions from the lowest band and we can 
drop the band index from our notation. The interband 
component is then 

^(k)=iV(l-|/ ,o(k)| 2 ). (29) 

Here we see that when the site overlap is negligible, the 
contribution from the higher bands is dependent only on 
the shape of the site Wannier function and provides no 
other information about the atoms in the lattice. The 
contribution to light scattering from the lowest band is 

S g (k) = |/ ,o(k)| 2 Y,(nAy k<Tj - T,) - (30) 

Deep in the Mott insulator regime the angular depen- 
dence of light scattering could be calculated using per- 
turbation theory, an approach that has been applied in 
theoretical work investigating Bragg scattering by Rey et 
al. @ and extended in Rist et al. However, perturb- 
ing J from zero leads to corrections to the unperturbed 
scattering pattern of order M(J/U) 2 and corrections to 
the unperturbed energies of order J @. When J < U 
these perturbations have negligible effects on the scat- 
tering pattern and the scattering pattern is well approx- 
imated by that produced when J = 0. We therefore 
proceed to calculate the scattering at finite temperature 
by assuming J = 0. 

To evaluate Eq. (|30|) in the Mott regime at tempera- 
ture T we need to calculate ^ u {^ u \h : jhi\'i)) u )e~ Eu / kBT / Z 
where \ip u ) are the eigenstatcs of the Bose-Hubbard 
Hamiltonian with energy E u . Making the approxima- 
tion J = [46j|, the eigenstates are simply number 

states |{n} u ) = \{n>j,j = 1,...,M}) with energies 

J2j Unj(rij — l)/2, where the total number of atoms, 

y"V rtj = Muq, is fixed with no a positive integer. At 
zero temperature the ground state is simply the number 
state with hq atoms at each site. 

For a translationally invariant system, the energy of 
these number states remains the same under any per- 
mutation of lattice sites. The eigenstates are then di- 
vided into degenerate groups labeled by v, where 
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FIG. 4. (Color online) Angular dependence of components 
of the static structure factor in the Mott insulator regime 
for a 30 x 30 lattice with no = 1 and V x = V v = V z = 
15Er. The temperature is 0m3E R /k B (U = QA2E R ) at 
which point the zero-temperature ground-state proportion is 
0.32. (a) Dependence on 9 of the structure factor components 
5 s o(k) (solid blue line), S 9 i(k) (dashed green line), and S&(k) 
(dotted red line) for <f> = 0. (b) As in panel (a) with cj> = 
7r/4. (c) Comparison of the Mott insulator structure factor 
(solid blue line) with that for the superfiuid (dashed green 
line). For comparison the superfiuid state is produced for the 
same lattice parameters as the Mott insulator, but with the 
s-wave scattering length scaled by a factor of 0.0006, giving 
a condensate population No = 807. 



is a representative eigenstate. The other members of the 
group are got by permutations P of the lattice sites, re- 
sulting in g v different states |P{n}„) each with energy 
E v . We can use this symmetry to evaluate the required 
matrix elements of Eq. ([5D)l . For each v we have 

{9v(n 2 ) v j = l 

g v (N 2 - M(n 2 ) v ) 
M(M - 1) 3 * ' 
(31) 

where (n 2 ) v = jj yj . (n>j ) is the average over the de- 
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w 



generate group v. This then gives 

3,1 P 



9v 



(W),, - Tin 



(M - 1) 



where 



i,i je{n:,y,z} 



(M 2 -^(k)) 



sin 2 (A/ J -k :( a/2) 
sin 2 (kj a/2) 



(32) 



(33) 



is the classical diffraction pattern from an M x x M y x A/ 2 
array of apertures with inter-site separation a in each 
dimension. 

Finally we find that the lowest-band contribution to 
the light scatter results in two terms. The first is the 
classical pattern due to the average density 



5 ff0 (k) = |/ ,o(k)| 2 n 2 ^(k), 



(34) 



and the second term is due to the number fluctuations 
on each site at finite temperature 



S 9 i(k) = 



|/ , (k)| 2 (M 2 -J-(k)) 



(M - 1)Z 



-EJk B T ( ,2 



((n 2 ),-n 2 ). (35) 



Neither of these terms result in excitation of the many- 
body state and the light scattering within the lowest band 
is purely clastic. 

For a number of sites of around 50 this formalism al- 
lows us to calculate the light scattering resulting from all 
possible excitations. We soon see however, that in the 
temperature range where the zero temperature ground 
state, \{n,j = no, j = 1 , • ■ • , M}}, is still significantly pop- 
ulated, the only states that play a significant role deter- 
mining the scattering distribution are the states involv- 
ing particle-hole excitations. These are the states that 
have no + 1 atoms at the sites pi,P2, ■ ■ ■ ,Pv and an equal 
number of sites hi, h%, . . . , h v with n<j — 1 atoms. There 
are g v = M\/[{M — 2v)\(v\) 2 } states with v particle-hole 
pairs and these have energy E v = vU and site number 



fluctuations of ((n 



2v/M. Restricting our 



calculations to just these states allows for calculations 
involving much larger lattices. 

In the limit T = the light scattering is purely due 
to the classical diffraction pattern and interband scat- 
ter, no excitation of the lattice occurs within the low- 
est band, a major difference from the superfluid phase. 
We can see why this is by considering the simple two 
site case. Switching from the site basis to the lattice 
momentum basis using 6 q = ^= j h e * q rj '> we see 
that the Mott state with filling factor one is Sj^lO) = 
l/2((^) 2 - (S^) 2 )|0), where q = f (r 2 - n). Light scat- 
tcring involving the wave-vector change k = q occurs by 



Stop 



FIG. 5. (Color online) (Left) Scheme for collecting photons 
scattered from atoms in an optical lattice to measure the tem- 
perature of the gas. A stop placed in the center of the Fourier 
plane of the first lens blocks the central diffraction peak. A 
second lens collects the unblocked photons which are then 
detected. (Right) Outline of the stop on diffraction pattern. 



two routes, b]J>o and b^bq. Each of these routes takes the 
Mott state to the same state, but with opposite signs, 
and the amplitudes cancel. The lack of excitations hence 
is a result of interference due to correlations in the lattice 
momentum distribution. 

In Fig. 2] we plot the components of the static struc- 
ture factor for a Mott insulator at finite temperature. 
The classical diffraction pattern remains the central fea- 
ture and does not change with temperature. The in- 
terband scattering leads to a background level of scat- 
tering around the classical diffraction peaks that scales 
with N, which may explain part of the background scat- 
tering observed in Ref. [TtJ that has the same scaling. 
The interband scattering is dominant over the scattering 
resulting from the site number fluctuations, which dis- 
appear as T — > 0. In Fig. [4jc) we compare the static 
structure factor for the Mott insulator with that for the 
superfluid. The stronger inelastic scattering in the super- 
fluid case leads to higher scattering intensity away from 
the diffraction maximum, although the difference is not 
as great as suggested in Lakomy et al. [H| due to the 
presence of the interband scattering component. 



IV. TEMPERATURE DEPENDENCE OF LIGHT 
SCATTERING 

In all the cases described above components of the 
static structure factor depend on temperature and we can 
use this dependence to observe changes in tem per ature 
[ljl[f6|]. In their recent paper, Ruostekoski et al. |l6| pro- 
posed using non-resonant light scattered from fermions 
in optical lattices to determine the temperature of the 
Fermi gas. It is important to test whether their method 
could also be applied to bosons in optical lattices, as an 
accurate in situ thermometer would be a useful experi- 
mental tool [13 • Furthermore Ruostekoski et al. neglect 
excitation of atoms into higher bands in their treatment 
of light scattering, and the effect of interband scattering 
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FIG. 6. (Color online) (a)-(c) Temperature dependence of the total number of photons scattered over all angles (solid line), 
the number of inelastically scattered photons (dot-dashed line), and the number of photons scattered due to the higher bands 
only (dotted line) for three many-body states on a 150 x 150 lattice, (a) Fermions in a lattice with V x = V y = 8Er, Vz = 15-Er 
and filling factor / = 0.5. (b) Superfluid state in a lattice with V x = V y = 3Er, V z = 20Er and filling factor / = 1. (c) Mott 
insulator state with V x = V y = V z = 20Er and filling factor / = 1. Note that for the Mott state the interband scattering makes 
up all the inelastic scatter, (d)-(f) Number of photons collected per second by the lens system described in the main text for 
the same lattices in panels (a)-(c). We also plot on the right axes, with a dashed green line, (d) Nfs/N, where Nfs is the 
number of atoms with energies below the Fermi energy; (e) No /N, where iVo is the condensate population for the superfluid; 
and (f) Nmi/N, where Nmi is the proportion of atoms in the zero temperature ground state for the Mott insulator. The 
scattering rates are determined for the D2 line of 40 K in the fermion case and for the D2 line of 87 Rb in the boson cases, each 
with a detuning of 2O7 from resonance and a laser intensity of U n = 5Wm~ 2 . 



on the temperature measurement must be quantified if 
this method is to be used in experiments. 

Following the scheme of Ruostekoski et at, scattered 
photons can be collected using a system of two lenses as 
shown in Fig. [5j where the number of photons detected 
gives a signal dependent on temperature. As noted by 
Ruostekoski et ah, the elastically scattered photons have 
little temperature dependence and placing a stop at the 
center of the Fourier plane of the first lens will improve 
the signal-to-noise ratio, by blocking the central diffrac- 
tion peak. 

In Figs.[BJa)-[njc) we plot the temperature dependence 
of the total number of photons scattered per second over 
all angles for our three many-body systems. Tempera- 
ture dependence is weak in all cases. We also plot the 
proportion of scattering that is inelastic, which is largest 
for the fermion case and smallest in the Mott insulator. 
In Figs.[6jd)-[6jf) we plot the temperature dependence of 
the number of photons collected, N C (T), by the lens sys- 
tem described above for the three many-body systems. 
We have assumed that the lens system has a numerical 
aperture sin# = 0.5 and the stop blocks light scattered 
with 6 < 0.06. In each case we see an increase in photon 
number with temperature that could be used to measure 
the temperature. 



As discussed by Ruostekoski et ai, inelastic scatter- 
ing processes heat the sample, and to measure the tem- 
perature without perturbing the system significantly we 
must limit the number of inelastic events IF in a sin- 
gle experimental realization. Following Ruostekoski et 
al. we take W = 0.1N. To measure the temperature to 
a useful accuracy it is then necessary to do multiple rep- 
etitions t of the experiment. The accuracy of the mea- 
surement then depends on the Poissonian fluctuations 
^/tN c {T) of the number of photons collected in the r 
repetitions. To achieve an uncertainty AT in the tem- 
perature measurement at temperature T we then require 
t(N c (T + AT) - N C {T)) ~ ^tN c (T) 0. 

For the experimental parameters discussed in Fig. [6j 
determining the temperature with accuracy of T = 
0.05T F (T F = 0.12E R /k B ) in the fermion case would 
require 39 repetitions at T = 0.5Tp and 7 repetitions 
at T = O.lTp. The presence of the inelastic interband 
scattering, which is largely temperature independent, re- 
duces the efficiency of the measurement, where for com- 
parison if we had neglected the interband scattering the 
required number of repetitions would have been 22 and 
4. In the superfluid case determining the temperature 
to accuracy Q.QQ2En/ks would require 23 repetitions at 
k B T = 0.02E R and 197 repetitions at k B T = 0.004T F . If 
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intcrband effects were excluded the number of repetitions 
required would have been 17 and 138, respectively. The 
temperature independent interband scattering causes the 
greatest reduction in efficiency in the Mott insulator case, 
where a temperature measurement at T = 0.027 'E^/ks 
would take around 26000 repetitions to achieve an ac- 
curacy of 10%. This results from the weak temperature 
dependence of the collected photon number relative to 
the amount of interband inelastic scattering when the 
zero temperature Mott ground state remains significantly 
populated, as shown in Fig.[5Jf). In contrast without the 
interband contribution all the remaining scattering in our 
model is elastic and in the limit U / J — > oo we would only 
need one repetition. 



The reduction in efficiency is due mostly to the in- 
creased heating resulting from the interband scattering 
events rather than to a decreased signal-to-noise ratio. 
This means that filtering out interband photons before 
collection (as may be possible due to their separation in 
energy by at least the band gap from the other scattered 
photons) does not increase the measurement efficiency 
greatly. In the Mott case above, if we filtered out the in- 
tcrband photons before photon detection the number of 
repetitions required would be reduced to around 24000. 
Similarly for the fermion and superfluid cases above the 
increase in efficiency is 5% or less. 



V. CONCLUSIONS 

We have shown how correlations between atoms in op- 
tical lattices are reflected in the way these systems scatter 
light. For the Fermi-Dirac gas we have seen that Pauli 
blocking leads to angular dependence of the inelastic light 
scattering, while for a bosonic superfluid a different an- 
gular dependence occurs as a result of interference and 
enhancement between scattering channels. In the bosonic 
Mott insulating state inelastic scattering is strongly sup- 
pressed by correlations between the atoms in momentum 
space. We have demonstrated that interband scattering, 
which has been neglected in other treatments, can sig- 
nificantly change the scattering patterns in a way that is 
largely independent of correlations within the lattice and 
temperature. 

We have used our theory to test whether light scatter- 
ing can be used thermometer for ultracold atoms, 
and have shown that intcrband scattering reduces the 
efficiency of temperature measurements based on light 
scattering, particularly in the Mott insulator case. Our 
results suggest that using light scattering to measure the 
temperature of the Mott insulator state will be ineffec- 
tive, while for the Fermi-Dirac gas and the superfluid this 
remains a feasible method. 
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